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We compare predictions of the Capillary Wave Model with Monte Carlo results for the energy gap and the 
interface energy of the 3D Ising model in the scaling region. Our study reveals that the finite size effects of these 
quantities are well described by the Capillary Wave Model, expanded to two-loop order (one order beyond the 
Gaussian approximation). 



An effective model widely used to describe 
a rough (fluid) interface is the Capillary Wave 
Model (CWM). In its simplest formulation one 
assumes an effective Hamiltonian proportional to 
the variation of the interface's area with respect 
to the classical solution. Until recently [1] the 
CWM has been studied only in its quadratic ap- 
proximation. The improvement of Monte Carlo 
(MC) calculations reached in the last years now 
allows tests of the CWM beyond the Gaussian 
level. This paper is a very short version of ref. 
[2]. We compare 2-loop predictions of the CWM 
for the tunneling mass gap and the interface en- 
ergy with MC results for a rough interface in the 
scaling region of the 3D Ising model. 

1. THE MODELS 

1.1. The Ising Model 

We consider the 3D Ising model on cubic lat- 
tices of size LijL'z in the x—, j/— directions, with 
L'z > Li, and of size t in z-direction. In the 
X—, J/— direction we always have periodic bound- 
ary conditions, whereas in z-direction either pe- 
riodic or antiperiodic boundary conditions are 
used. The Hamiltonian is a sum over nearest 
~Y^<ii> ^i *i' ^nd the 



neighbor pairs, H 



Boltzmannian is exp(— /3_ff). While in inflnite vol- 
ume, for (3 > (3c, the system shows a spontaneous 



symmetry breaking, in flnite volume this cannot 
occur, and interfaces spontaneously appear, sep- 
arating domains of different magnetization. The 
critical coupling of the model is (3^ = 0.221652(3) 
[3]. In the region between the critical and the 
roughening coupling (3^ = 0.4074(3) [4] the inter- 
face behaves essentially as a 2D critical system. 
It is a widely anticipated assumption that this 
system can be described by a (continuum) model 
with a Hamiltonian that is proportional to the 
area of the interface. 

1.2. The Capillary Wave Model 

If the interface is described by a single valued 
height function 4> over the torus with coordinates 
{) < Xi < Li, i = 1,2, the CWM Hamiltonian is 
'H(<j>)=A(<j>)-LiL2,mth 



A\ 




(1) 



We are interested in the free energy F, deflned 
through 



F = - keT InZ 

Z = X exp{-aLiL2)Zq{a,Li,L2) 



Z„ 



{D<^\ exp {-n 



(2) 
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A is an undetermined constant in this approach. 
By expanding the root in eq. (1), one can set up 



a perturbative expansion for Z^. The expansion 
parameter is {(jA)~^, with A = L\ L'j. Unfortu- 
nately, the model is UV (by naive power counting) 
nonrenormalizable at large momenta. However, 
at least to 2-loop, the same finite results are ob- 
tained by a large class of regularizations [5,2,6]. 
To 2-loop order, one gets Zq = Zq Zq , with 



1 



Z(")(„)= ^(i„)/^(i) 



(3) 



r\ is the Dedekind eta function, and u = L^j L\ is 
the asymmetry parameter. For the 2-loop contri- 
bution one obtains 



Zf) = ! + /(«) (^A)-i 



i 
with 



(4) 



1 \\T^ 



/ («) = 2 I [ 6 "^2 {iu)\ - -uE, iiu) + 4 I . (5) 

Here, E2 denotes the first Eisenstein series. Cor- 
rections from higher loops are of order (aA)~^. 

2. OBSERVABLES AND MC 

The observables that we used to check predic- 
tions of the CWM were the energy gap AE and 
the interface energy Es- 

Due to tunneling, there is a small energy gap 
AE between the symmetric and antisymmetric 
ground state of the Ising model in a finite geom- 
etry. In the dilute gas approximation, the energy 
splitting is directly proportional to the interface 
partition function, AE = Z(a, Li, L2), and thus 
allows for a comparison with the CWM partition 
function. We computed the energy gap for sev- 
eral values of /3 and for various combinations of 
Li and L2 employing two different methods: the 
time slice correlations method (TSC) (see, e.g. 
[7]) and the boundary fiip method (BF) [8]. 

Our second observable was the interface energy, 
that can be computed from the CWM partition 
function, Es = -Z'^ {dZ/dp). The 2-loop result 



Es{a, Li,L2) = a'LiL2 - ^ - ^S/jZf') . (6) 



A ^m P " 



For the Ising model, we estimated the interface 
energy from the difference Es = {H)p — {H)a of 
the energy expectation value with periodic and 
antiperiodic boundary conditions, respectively. 
These expectation values were computed using a 
microcanonical demon algorithm in combination 
with a particularly efficient canonical update of 
the demons [9]. The algorithm was implemented 
using the multi-spin coding technique. 

3. MC DATA ANALYSIS 

3.1. Energy Gap 

Our /3- values are listed in table 3.1. For these 
values, the infinite volume bulk correlation length 
ranges from ^ Ri 2.6 for j3 = 0.2275 to ^ Ri 4.5 
for P = 0.2240. We simulated for typically 10 
to 20 different combinations of Li and L2, where 
L2 > Li. We always kept Li larger than the 
inverse deconfinement temperature of the dual 
gauge model, t was usually fixed to 120, for sim- 
ulations with very large Li and L2 also larger. 
To test the CWM prediction, we fitted the 2-loop 
result to our MC data. Note that there are only 
the two parameters A and a to be fitted. 



Table 1 

(3-values for AE -measurements. T means TSC 

method, B means BF method 



P 



T mm 



0.2240 18 
0.2246 13 



B 

T 



P 



T mm 



0.2258 
0.2275 



12 
10 



T 
T,B 



Here, a' = daa and A' = daX 



^P 
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As an example, we show in fig. 1 our MC results 
for three of the /3-values, together with the best fit 
curves. Obviously, the size and shape dependence 
of AE is nicely described by the 2-parameter fit 
with the 2-loop CWM partition function. 

Compared to fits with the classical result, the 
reduced x^ of the fits with the 2-loop-result were 
always more than one order of magnitude smaller. 

3.2. Interface Energy 

For the interface energy, the 2-loop contribu- 
tion appears as an additive correction to the 



Gaussian (and classical) result. We made our cal- 
culations at /3 = 0.24. Precise estimates of a and 
it' were known from other MC experiments. We 
considered these estimates as external input. Let 
us define the following energy differences: 



AEs (L) = Es (2L, L/2) - Eg (L, L) 



(7) 




Fig 1: MC data and best 2-loop fits for dif- 
ferent values of (3 versus A = L1L2. From 
top to bottom, the first three lines correspond to 
j3 = 0.2246, the next three to j3 = 0.2258, and the 
last two to j3 = 0.2275, 



The CWM prediction for this quantity is 



A^' 



(CWM) 



(721.2 



[/(4)-/(l)] 



(8) 



In fig. 2, we plot this prediction together with the 
corresponding estimates obtained from our MC 
data for the interface energy. A very nice agree- 
ment is seen for _L > 16. 

CONCLUSION 

Our results show that the finite size effects of 
interface properties in the scaling region of the 
Ising model are well described by the 2-loop ex- 
panded CWM. 




Fig 2: Comparison between MC and CWM pre- 
dictions for AE's, The dotted line is the theoret- 
ical prediction if 2-loop corrections are neglected. 
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